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Abstract 

We consider an application of pooled stepped chutes where the 
transport in each pooled step is described by the shallow-water equa- 
tions. Such systems can be found for example at large dams in order 
to release overflowing water. We analyze the mathematical conditions 
coupling the flows between different chutes taken from the engineering 
literature. We present the solution to a Riemann problem in the large 
and also a well-posedness result for the coupled problem. We finally 
report on some numerical experiments. 
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1 Introduction 

This work deals with water behavior of so-called pooled stepped chutes. 
This geometry frequently appears in the real water dams and it also appears 
in mountain rivers to control the bed load transport. In both cases the 
main concern is to spill excessive noodwater in additional channels next to 
the dam structure. These are called pooled stepped chutes or pooled steps. 
Within the pooled steps additional weirs perpendicular to the flow direction 
are introduced to increase energy dissipation. This problem has been gained 
some attention in recent years in the engineering community, see e.g. [UE1E1 
IHIIUII2 H 124} 126]. However, only a few mathematical discussions are currently 
available [16} I23j . In particular, the modeling and design of the spillways 
and stepped channels have so far been addressed using experiments and data 
fitting techniques, see e.g. |26j . From the measurements empirical formulas 
have been derived and used in sophisticated simulations. The measurements 
taking into account complex geometries as well as material properties and 
the air-water mixture leading to a variety of different empirical formulas 
and tables [TTJ. 
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So far, the mathematical discussion has been limited to a consideration 
of the effect of the weir at the end of a spillway neglecting the dynamics 
of the water inside the pooled channels. Here, we discuss the mathematical 
implications of considering the coupled problem, i.e., the dynamics inside 
the pooled steps and the (empirical or theoretical) conditions imposed closed 
to the weir. Typically, the water flow in the channels is described by the 
shallow-water equations whereas the effect of the weir is given by some al- 
gebraic condition. We treat this problem using a network approach with the 
conditions at the weir as coupling conditions. We present a well-posedness 
result for a simple condition based on energy dissipation. The recent liter- 
ature offers several results on the modeling of systems governed by conser- 
vation laws on networks. For instance, in [TJ [2] [TO] [T3] the modeling of a 
network of gas pipelines is considered. The basic model is the p-system or, 
in |14[ 118] . the full set of Euler equations. The key problem in these papers 
is the description of the evolution of fluid at a junction between two or more 
pipes. A different physical problem, leading to a similar analytical frame- 
work, is that of the flow of water in open channels, considered for example 
in [2D]. 

Consider a water flow in an open canal affected by a weir or small dam 
at the point x = 0. If the water level becomes greater than the height of the 
weir, some water passes over the weir. Similarly to the models in [15 t I19 |. [20] 
we describe the dynamics of the water by the shallow water equations, while 
the interaction with the weir is described by coupling conditions. 

Let (h,v)(t,x) be respectively the water level (with respect to the flat 
bottom) and its velocity for i / 0. The shallow-water equations in each 
canal are given by 

d t h + d x (hv) = 
' dt (hv) + d x (hv 2 + \gh 2 ) =0 X * °' (L1) 

where g is the gravity constant. Two canals are coupled by so-called pooled 
steps [26j . A sketch of this situation is given in Figure [TJ In the engineering 
literature [9] [25] [26] (see also Remark II. ip water of height h~~ flowing over 
a weir of height H~ generates a flow Q 

( \ 3/2 

Q = C Jg-(h- -H-) , (1.2) 

where C is a constant depending on the air-water-ratio and the detailed 
geometry. In many situations C = 0.6 is used. The equation (|1.2p is called 
3/2— law. Formally, it can be derived from the following idea: the potential 
energy of the water flowing over the weir is transformed to kinetic energy. 
Hence, we have the balance 
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Figure 1: This figure represents two connected pooled steps with a weir in 
between and the indication of the various heights 



Solving for v in terms of h we obtain the previous formula and 1 — C is the 
percentage of energy loss during the change of potential to kinetic energy. 
We use this idea to deduce the coupling condition for (jl.ip at x = 0. We con- 
serve the total water over the weir and hence h + v + = h~v~ . The difference 
in the amount of water overflowing the weir is \h~ — H~] + — [h + — H + ^ + . 
This defines the velocity (and its sign) at the weir according to the balance 
of potential and kinetic energy. Hence, the coupling conditions are 

W = C ([h--H-] + - [h+ - H+] + ) ■ J[h--H-} + -[h+-H+} + 

< 

h + v + = h~v~, 

(1.3) 

where C = 0.6^, (h~ ,v~) = (h,v)(t,0-), {h + ,v + ) = (h,v)(t,0+) while 
are the heights of the weir to the left and the right (see Figure [T]) . 

Remark 1.1 Obviously, hi. 2^1 is only a first approximation on the complex 
dynamics at the weir, see |7| \23\ \26\/ . As outlined in the introduction there 
exists a variety of empirical formulas in the engineering community. Many 
of them include further effects as for example the water-air ratio of the over- 
spill or the roughness of the channel bottom. For example in [26, Equation 
7. 7] the following relation has been determined 

/ \3/2/2 , — 

Q=(h--H-j (-//v/25 

/ \3/2 / n5/2 

C 1 ^/g(h--H-) +C 2 Vg{h--H-) /H~. 

Here, we have /i = 0.611 + 0.08(/i~ — H~)/H~ and the constants are C\ = 
| \/2 0.611 and C2 = | V% 0.08. Since the coefficient C\ is roughly eight 
times larger than C2, this equation is very similar to hl.2\) . Another example 
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is given in J3j/ ; }26\. Equation 2.50-2.51], where the following formula has 
been proposed for the flow with H = h~ — H~ 



(v-) 2 

Q = v~H(k c + k d ) = 0.15 - 0.45 + 

2gH 

( (v-f V ( ({v-f 

0.57 - 2 V-rr - 0.21 exp 10 - 0.21 

\2gH / I 2gH 




V 

The case of no weir, i.e., H~ = and H + > 0, lead to the following 
studied formulas, e.g. f23jj, f2b\ Equation 2.38] 



Q = TT0.715/T 
or JEEi, JM Equation 2.39-242] 

, / , \ 1.275 

Q = V -H + [h-/H + \ 

We restrict our discussion to the still commonly used hl.2\) to outline the 
ideas. Further note that additional empirical formulas for the arising wave in 
the outgoing pooled step are not needed in our approach, since these dynamics 
are fully covered by the shallow-water equation. 

2 The Riemann problem for a single weir 

By Riemann Problem at the weir we define the problem (jl.ip . (jl.3p with 
initial data 

<M) = '""I (2.4) 

I {n r ,v r j tor x > 0. 

Definition 2.1 A solution to the Riemann Problem (jl.ip . (|1.3p . (|2.4p is a 

function (h,v) : M + xl->l 2 such that (t,x) — > (h,v)(t,x) is self-similar 
and coincides in x > with the restriction of the Lax solution to the standard 
Riemann Problem for (|1.1|) with initial data 

(M) = {(M)(*,0+) for,< 
I (h r ,v r ) for x > 0. 

«;/n/e coincides in x < uiii/i i/ie restriction to of the Lax solution to the 
standard Riemann Problem for (jl.ip with initial data 

I (/i,v)(i,0-) for x > 0. 



moreover 



(h,v)(t,0±) = (/i ± ,w ± ) sate/y (OP . 
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For studying the Riemann problem we first collect the standard expres- 
sions for the eigenvalues, eigenvectors and Lax curves for the shallow water 
equations (|l.ip . 

The 2x2 system of conservation laws in (jl.ip has the eigenvalues Ai, A2 
and the eigenvectors r%, r 2 , where 



Ai(/i, v) = v 

n{h,v) = 



X 2 (h,v) = v + y/gh 



1 



-v + y/gh 



r 2 (h,v) 



1 



v + y/gh 



(2.7) 



V( ftift „)Ai(/t,u) • r i(M)^yf > V (hM) X 2 (h,v) ■ ^(/i,^)^^! > 0. 
The Lax curves of the first and second family, described in Figure [21 are: 



v = £~t(h;h ,v ) 



v = C^(h;h ,v ) 



v -2 (y/gh - y/gh^) h < h 

vo-(h-ho)y/y^ h>h , 

v + 2 (y/gh - y/gh^) h > h 

vo + (h-ho)y/y^ h<h . 



(2i 



(2.9) 



The reversed Lax curves of the first and second family are given by: 



v = C 1 (h;h ,v ) = < 



v = C 2 (h;h ,v ) 



v - 2 (y/gh - y/gh^) h > h 



v -(h-h )^Jy^ h<h , 

v + 2 (y/gh- y/gho) h < h 



(2.10) 



(2.11) 



Concerning the coupling condition, we study the states which satisfy it at 
the point x = 0. We restrict ourselves to the subsonic region \v\ < y/gh. 
From (|1.3p . since the heights are always positive, it follows 



Moreover 



signw 



sign 



C 



H 



}/ 



H 



+ 



which becomes 

I h ~ \ v ~ 



h+ 



sign v , 
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Figure 2: Lax curves for the shallow water system in the (h, hv) plane 



h + - H 



h~ - H' 



h \v 
C~ 



• sign v 



We have also to add the equality of the fluxes, hence we obtain 



[h+-H+] + = [hr-H- 
h + v + = h~v~ . 



h-\v~\ \ 3 

^r— 1 ] • sign V 



(2.12) 



Consider first the case where the water level to the right is below the weir, 
h + < H + so that \h + — H + ~\ + = 0. If also h~ < H~ then no water crosses 
the weir, while if h~ > H~ the water crosses the weir flowing from left to 
right. The left state must satisfy v~ > and 



JET 



h \v 
C~ 



o.. 



which can be rewritten as 



h v 



C h 



(2.13) 



The curve {h~ ,v~) where h~ , v~ satisfy ()2.13p consists of all the left states 
which can be connected to a right state with h + < H + (see Figure [3|). 
We call the support of this curve T u . The case h + > H + , h~ < H~ is 
symmetric, hence we call Ti the support of the curve 



-h+i 



C h + -H 
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Since C [h~ - H~] \ = CJg [h~ - H~] \ < y/g (hr) ° both T u and F t lye 
in the subsonic region. If both levels are above the weir, then (|2.12p gives 
a unique state (h + ,v + ) which connects a given state (h~,v~). We call this 
function (h + ,h + v + ) = &(h~ , h~v~). We have the following lemma. 

Lemma 2.2 Consider the (h, hv) plane, and define the following sets de- 
picted in Figured 





( 

= Uh,hv) 




= {(h,hv) 


A~ 


= 


B~~ 


= {(h,hv) 


n, 


= |(/i,H 




= {(h,hv) 


A+ 


={(h,hv) 


B + 


= l(h,hv) 



h-H~ 



h-H~ 



h-H' 



-C 



h-H' 



< hv < 



Then, we have that the coupling condition induces the following assertions: 

• no left state (hr ,h~v~) G $7 M can be connected to any right state 
(h + ,h + v + ) and no right state (h + ,h + v + ) G Qi can be connected to 
any left state (h~ , h~v~); 

• any left state in T u can be connected to any right state in S n with the 
same flux: h~v~ = h + v + ; 

• any right state in T[ can be connected to any left state in Y>i with the 
same flux: h~v~ = h + v + ; 

• any left state in A" can be connected to one and only one right state 
in A + and any right state in A + can be connected to one and only one 
left state in A" ; 

• any left state in B~ can be connected to one and only one right state 
in B + and any right state in B + can be connected to one and only one 
left state in B~ . 
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Figure 3: Representation of the regions introduced in Lemma 12.21 

We omit the proof since it is a straightforward analysis of the coupling 
conditions (|27L2l) . 

We call : A~ U B~ — > A + U B + the one to one map that associates to 
any left state in A~UB~ the corresponding right state in A + U5 + connected 
through the coupling condition: 



in such a way that the left state {h + ,h + v + ) = ^(h~,h' ~v ) is connected 
through the coupling condition to the right state (h~ ,h~v~) for any 
(h~,h~v~) € S M U A~ U B~ . 

We are now able to show that the Riemann problem can be solved in 
the large. 

Proposition 2.3 For all states (h[,vi) and (h r ,v r ) in the subsonic region, 
the Riemann problem (|1.3|) . (|2.4|) has a unique solution satisfying 

Definition \2.1l whose constant states depend continuously on the initial data 
and are constructed glueing together a wave of the first family, the coupling 
condition and a wave of the second family. 

Proof. Given two states (hi,vi) and (h r ,v r ) in the subsonic region, we 
proceed in the following way for solving the Riemann problem. Draw the 
Lax curve of the first family in the (h,hv) plane through (hi,vi). Since, in 
the subsonic region, it is strictly decreasing, it intersects in one and only 
one point (with h > 0) the convex and non decreasing curve T u . We call the 




We also define 




for any {h, hv) € £/ 
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unique point of intersection (h*,v*). Then, we consider the non increasing 
curve 

(^h,h*v*) foxh<h* 
j(h) = < V / ; x (2.14) 

^(h,hjOf(h;hi,vi)) for h > h*. 

Next, we take the inverse Lax curve related to the second characteristic field 
passing through the right state (h r ,v r ). This curve is strictly increasing in 
the subsonic region and in the upper supersonic region, therefore it has one 
and only one point of intersection with 7 (for positive h), denoted by (k, w). 
This point might belong to the upper supersonic region. The Riemann 
problem is finally solved in the following way: take the subsonic (or lower 
supersonic) state (k*,w*) on the first Lax curve such that kw = k*w* and 
connect (hi,vi) to (k*,w*) with a wave of the first family. This wave has 
negative velocity since the curve belongs entirely to the subsonic region or to 
the subsonic region and the lower supersonic one. Then, the points (k* , w*) 
and (k,w) satisfy the coupling conditions (|l,3p . Finally (k,w) and (h r ,v r ) 
are connected by a wave of the second family which travels with positive 
velocity even if (k, w) happens to be in the upper supersonic region, see 
Figure H □ 



3 A well posedness result 

This section is devoted to show the well-posedness of the Cauchy problem 
for system (jl.ip for initial data around a constant subsonic state satisfying 
the coupling condition (|1.3|) . For simplicity we suppose that the water level 
to the right of the weir is below the weir level, that is h + < H + . This implies 
that the water can only flow from the left to the right of the weir when its 
left level overflows the weir. 

Theorem 3.1 Given Hq , H + and two constant states in the subsonic re- 
gion (ho,vo) and (hi,vi) such that 



h v = C 



Ho 



hoVo = h\Vi, ho > H , h%< H + , (3.15) 



then, there exists a closed domain 



T>cUh,v)e (/ l o^o)x(-oo,o) + (^i^i)x(o,oo ) + (L 1 nBV) (m ; m 2 ) j 

containing all functions with sufficiently small total variation in x > and 
x < and semigroups 



defined for all H sufficiently close to H , such that 
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Figure 4: These figures describe how to solve two Riemann Problem (first 
and second row). The states on the left and on the right of the weir are rep- 
resented respectively on the left and right figures. The green lines separate 
the subsonic and supersonic regions 



1 ) for all t, s > and u 6 T> 

Sq u = u, 

2) for all u,v S T>, , in a suitable neighborhood of Hq and t,t' > 0: 



Sf 1 u - sf 2 v 



< L- \\\u- w|| L i + \t-t'\ + t ■ \H{ - Hz\\ ; 



3) if u £ V is piecewise constant, then for t small, S^ 1 u is the glueing of 
solutions to Riemann problems at the points of jump in u and at the weir 
in x = 0; 



4) for all u G V, the map u(t,x) = \ Sf u J (x) is a weak entropy solution 
to CO]), (USD (see Definition 4-lJ, Jjj\ Definition 2.1]). 

S is uniquely characterized byUty, 0) and\3$. 
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Proof. Following [121 Proposition 4.2], the 2x2 system (jl.ip defined for 
x £ K can be rewritten as the following 4x4 system defined for x G M + : 



d t U + J r (U) = (t,x)€ 
b(U(t,0+)) =g(t) t£R+. 



(3.16) 



the relation between U and u = (h,v), between T and the flow in (jl.ip 
being: 



U(t,x) 



h(t, —x) 
-(hv)(t, —x) 

h(t, x) 
(hv)(t, x) 



T(t,x) 



U 2 

% + but 

Ua 



with x £ M + ; whereas the boundary conditions becomes 



(3.17) 



g(t> 



H- - H Q 




b(U) 



Ui 



2 

U 2 \ 3 



c 

Ua + U 2 



The thesis now follows from [11[ Theorem 2.2]. Indeed the assumptions 
(7), (b) and (f) are therein satisfied. More precisely the eigenvalues of the 
Jacobian of the flow in (|3.17[) are 

Ui Us 

therefore in the subsonic region exactly two are positive and exactly two are 
negative. Since here ^(t) = 0, j(t) = 0, condition (7) in [ITJ Theorem 2.2] 
is satisfied with 1 = 2. Concerning condition (b), the positive eigenvalues 
with the corresponding eigenvectors evaluated at U = (ho,—hoVo,hi,hiVi) 
are 

/ 1 \ 



-Vo + y/gho, R3 



-vo + y/gho 





A4 = v\ + y gh 



1; 



R \ 






1 



\ 



\v 1 + yfgKl J 



Conditions (|3.15p imply b ( U ) =0, and 

Db(U)R 3 , Db{U)RA 



-v + y/gho 





vi + \[gh\ 
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The determinant of the above matrix is given by 




«0 • -T7 



V\ + 



which is strictly positive since (/io,fo) an d (h\,vi) belong to the subsonic 
region. Thus condition (b) is satisfied. Concerning condition (f), system 
(|3.16p is not necessarily strictly hyperbolic, for it is obtained glueing two 
copies of system (jl.ip . Nevertheless, the two systems are coupled only 
through the boundary conditions, hence the whole wave front tracking pro- 
cedure in the proof of Theorem 2.2] applies. Concerning the Lipschitz 
dependence on the height H~ , observe that as in [11[ Theorem 2.2] we have 
for 



g(t)=H^-H { 



o > 



g(t) = H 2 -H i 



o > 



and hence 



/ \g{T)-g{r)\ dr = t 
Jo 



Hi — H 



□ 



4 Computational results 

We present some numerical results on pooled steps using a finite-volume 
method in the conservative variables to solve for the system dynamics in 
each canal. The coupling conditions (jl.3p induce boundary conditions for 
each canal at each time-step. For given data U := (h , (hv)^) close to the 
weir the conditions yield the boundary states at each connected canal. In 
the numerical computation of the boundary states we proceed as in Section 
[2] using Newton's method applied to (|1.3p where and (hv)^ are given 
by the forward (backwards) l-(2) Lax-wave curves through the initial state 

According to the solution of the Riemann Problem described in Section 
[2] there are three possible scenarios. If < the coupling condition (|1.3p 
reduce to v = and no water passes the weir. If h~ > H~ and h + < H + 
the water is overflowing the weir from the left to the right. If h + > H + 
and h~ < H~ the water is overflowing the weir from the right to the left. If 
h + > H + and h~ > H~ the water is flowing from the left or from the right 
depending on the sign of (h~ — H~) — (h + — H + ). 

Now, we present a numerical result for overflowing three connected pooled- 
steps. Each canal has length L = 1 and the height of each weir is 1.5. We 
simulate two situations. In the first case the initial water level in each canal 
is low (equal to H~) in the second case the canals are already full (height is 
equal to H + ). In both case the water initial is still v = and a wave with 
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a height h = 2 H + and hv = 5 enters on the first canal. This wave lasts 
until T = 20 and is then followed by a wave of height h = H~ and zero flux. 
The simulation time is T = 60 for both scenarios. We present snapshots 
of height and velocity at different times. The solid lines are the weir, the 
dotted line is ground level. It is a pooled step and therefore the first canal 
is on a higher level above ground than the last one. 

In the first scenario (see Figure E]) the wave enters the system of pooled 
steps and overflows the connected canals. After this wave passed the system 
slowly reaches again an almost steady state. In the second scenario (see 
FigureEJ) we observe initial dynamics due to the overflow of the pooled steps. 
On the second canal these dynamics interact with the incoming wave. After 
the wave passed the system slowly reaches a steady state with heights below 
critical (i.e., h < H + ) in the first and second canal. The water is still flowing 
over the third weir in this simulation. 
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University where part of this work was completed. 



13 



Figure 5: Simulation results for a strong wave entering the pooled step. 
Initially, the water level on each step is below critical. The solution (h,v) 
at different times is depicted. The time increases from left to right and top 
to bottom. 
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Figure 6: Simulation results for a strong wave entering the pooled step. 
Initially, the water level on each step is above the critical level. The solution 
(h, v) at different times is depicted. The time increases from left to right 
and top to bottom. 
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